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ABSTRACT 

This is a computational and experimental study for soot formation 
and radiative beat transfer in jet diffusion flames under normal gravity 
(l-g) and microgravity (0-g) conditions. Instantaneous soot volume 
fraction maps an: measured using a full-field imaging absorption 
technique developed by the authors. A compact, self-contained drop rig 
is used for microgravity experiments in the 2.2-sccond drop tower 
facility at NASA Lewis Research Center. On modeling, we have 
coupled flame structure and soot formation models with detailed 
radiation transfer calculations. Favre-averagcd boundary layer equations 
with a k-e-g turbulence model are used to predict the flow field, and a 
conserved scalar approach with an assumed /1-pdf are used to predict 
gaseous species mole fraction. Scalar transport equations are used to 
describe soot volume fraction and number density distributions, with 
formation and oxidation terms modeled by one-step rate equations and 
thermophoretic effects included. An energy equation is included to 
couple flame structure and radiation analyses through iterations, 
neglecting turbulence -radiation interactions. The YIX solution for a 
finite cylindrical enclosure is used for radiative heat transfer calculations. 
The spectral absorption coefficient for soot aggregates is calculated from 
the Rayleigh solution using complex refractive index data from a Drude- 
Lcmmtz model. The exponential-wide-band model is used to calculate 
the spectral absorption coefficient for H 2 0 and CO v It is shown that, 

when compared to results from true spectral integration, the Rosseland 
mean absorption coefficient can provide reasonably accurate predictions 
for the type of flames studied. The soot formation model proposed by 
Moss, Syed, and Stewart seems to produce better fits to experimental 
data and more physically sound than the simpler model by Khan et al. 
Predicted soot volume fraction and temperature results agree well with 
published data for a normal gravity co-flow laminar flames and turbulent 
jet flames. Predicted soot volume fraction results also agree with our 

data for l-g and 0-g laminar jet flames as well as l-g turbulent jet 
flames. 
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# Presently Design Engineer, Caterpillar Co.. Aurora, IL 60507. 


NOMENCLATURE 

b constant in the model for modified equilibrium temperature 
C ( ) rate constant for nucieation (a, 8), growth (/?), coagulation (y), or 
oxidation (x ) 

C* rate constant in the model by Khan et al., Eq. (9), C* = 16 8 
kg/Nms 

E s activation energy in the model by Khan et al., Eq. (9), E t = 

40000 cal/mole 
/ gas mixture fraction 

/ v . f v soot particle volume fraction, and f v - pjjp 
g variance off 

H total enthalpy 

k turbulent kinetic energy 

^ a length scale in the heat sink model 

m ( j exponent on the mole fraction in rate equations for nucieation (a) 
or coagulation (y) 

N, N soot particle number density, and N= N/(pn Q ) 
riQ Avogadro’s number, ^ = 6 x 10 26 

P f u partial pressure of unbumed fuel 

q r , q r radiative heat flux and radiative heat flux vector, respectively 
r radial distance from the flame axis 

S$ source term for property <p in Eq. ( l) 

T flame temperature 

P e q equilibrium (adiabatic) flame temperature 

P eq.max tnaximum equilibrium (adiabatic) flame temperature 
surroundings temperature 

T ( } activation temperature in rate equations for nucieation (a), 
coagulation ( 7 ), or oxidation (x ) 

u, v axial and radial gas velocity component, respectively 

v, radial direction thermophoretic velocity for soot particles 
axial distance from the nozzle exit 

X gas species mole fraction 

cc soot particle nucieation rate, in the particle transport equation for 
number density 

P soot particle surface growth rate 
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Y soot panicle coagulation rate 

S soot panicle nucleation rate, in the particle transport equation for 
volume fraction 

£ emissivity of soot/gas mixture 
e turbulent dissipation energy 

X soot particle oxidation rate 

P'jjj effective viscosity in Eq. (1) 

p mass density of gas mixture 

p s mass density of soot particles 

c Stefan-Boltzmann constant, a = 5.670 x 10“ 8 W/m 2 K 4 
<j> a general variable representing 1, u, f,k,e,g, / v , A\ or H 
< p local fuel/air equivalence ratio 

INTRODUCTION 

The subject of soot formation and radiation heat transfer in micro- 
gravity jet diffusion flames is important not only for the understanding 
of fundamental transport processes involved but also for providing 
findings relevant to spacecraft fire safety and soot emissions and radiant 
heat loads of combustors used in air-breathing propulsion systems. The 
longer flame residence time and the dominance of effects such as 
thermophoresis and diffusion, both as a result of eliminating buoyancy- 
induced convective flow, make microgravity an idea] environment for 
studying diffusion flame fundamentals. Comparisons between normal 
gravity (1-g) and microgravity (0-g) results, and between measurements 
and model predictions, will be used to improve our understanding about 
diffusion flames. Considerable amount of data have been published on 
microgravity jet diffusion flames (Bahadori and EdeLman, 1993). Most 
of these data are from photographic measurements, with some wide- 
angle radiometry, far-field temperature and gas species measurements. 
Our objective is to provide local measurements and modeling for soot 
volume fraction, flame temperature, and radiative heat fluxes in 
microgravity jet diffusion flames. 

In terms of our fundamental understandings of normal gravity 
diffusion flames, it is believed that gas-phase processes in laminar and 
turbulent flames are similar and relatively well known, with turbulent 
flames usually approximated as wrinkled laminar flames, known as the 
laminar flamelet concept (Bilger. 1976, 1977; Faeth and Samuelson, 
1986; Sivathanu and Faeth, 1990b). However, this similarity becomes 
somewhat conditional or invalid for soot and radiation properties (Gore 
and Faeth, 1986; Sivathanu and Faeth, 1990a; Sivathanu et al., 1990). 
Both soot formation and radiation heat transfer are inherently difficult 
subjects. Soot inception is much slower than gas-phase reactions, soot 
particles have much smaller diffusivities, and their growth, coagulation, 
and oxidation have not been well understood. The difficulties in treating 
detailed radiation heat transfer are to find accurate, yet computationally 
efficient, solutions for nonhomogeneous and nongray media, as well as 
accurate models for soot and gas radiative properties and for turbulence- 
radiation interactions. 

Our approach is to use simple, yet established, models for flame 
structure, which include Favre-averaged bound ary-layer-type governing 
equations with a k-c-g turbulence model for flow field and a conserved 
scalar approach with an assumed /3-pdf (probability density function) for 
gaseous species concentration profiles. Soot transport is described as 
that of a scalar property, including thennophoretic effects. Two existing 


soot formation models, with some improvements, will be tested and 
compared. Both models use one-step reactions to describe various 
mechanisms in soot formation. Radiative heat flux and its divergence 
are calculated using the YIX solution (Tan et al., 1990; Hsu et al., 1993), 
which is a numerical solution for the integral formulation of the radiative 
transfer equation. Total radiative properties, such as the Rosseland and 
the Planck means, are introduced to reduce computational times. Results 
based on these means are compared against results from detailed spectral 
integration. The flame structure solver and the radiation heat transfer 
solver are coupled through an energy equation, and they are solved 
iteratively until a set convergence is achieved. 

Predicted soot volume fraction and flame temperature results are 
compared to published data for 1-g co-flow laminar flames, as well as 
for 1-g and 0-g jet laminar and turbulent flames. 


MODELING OF JET DIFFUSION FLAME STRUCTURE 


The structure of turbulent jet diffusion flames is modeled using the 
Favre-averaged equations for conservation of mass, momentum, and 
mixture fraction. A conserved scalar approach (Bilger, 1976; Jeng and 
Faeth, 1984; Gore and Faeth, 1984) with an assumed probability density 
function (pdf) and a k-c-g turbulence model (Lockwood and Naguib, 
1975) are used. All governing equations can be written in a general form 
as (Faeth et al., 1985) 




1 d 


T W) * -j- r irpvf) = - r J- r [rp eff ., - + S t 


3<p 


( 1 ) 


where <p = 1, u (velocity), f (mixture fraction), k (kinetic energy), e 
(dissipation), or^ (variance of/). Details for fi e ^ ^ (effective viscosity), 
S$ (source) and assumptions are given in those references. Buoyancy 
effect is considered in mean flow only, neglecting buoyancy-turbulence 
interactions. For laminar flames, Eq. (1) is simplified accordingly. 

The system represented by Eq. (1) is solved using the block- 
tridiagonal code of Chen et al. (1987). State relationships constructed 
from equilibrium calculations using STANJAN (Reynolds, 1986) are 
introduced to decouple chemical reactions from flow calculations. More 
accurate laminar flamelet approaches (Bilger, 1977; Rogg, 1993) wiH be 
considered in the future. 


MODELING OF SOOT FORMATION AND OXIDATION 

We have considered two sets of soot formation and oxidation 
models. Both describe the transport of soot particles as that of a scalar 
property using Eq. (1). For soot transport in laminar flames, a radial 
thennophoretic velocity term, v f = - 054 (v / T)(JT / Jr) , is added. The 

two-equation model (Moss et al., 1988; Syed et al.,1990) is based on 
number density (N) and volume fraction (f v ), and the two respective 
source terms are 

SN=a-Pp 1 N 2 -r^ n xpf v - ]n N in , (2 ) 

P" o 

s /,= n o 1 (f-X)Pf v 2n N ln + 5, p = f v = BiA ■ (3) 

26 . P 
where Hq = 6x10 is Avogadro's number and p s - 1.8 g/cm 3 is mass 

density of soot. Soot oxidation models of Leung et al. (1991) and 
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Fairweather et al. (1992) are modified to include oxidation by both 0 2 
and OH. Rate -equations of Arrhenius type are used to model nucleation 
(a), growth (0), coagulation (# and oxidation (*). These rate equations 


are given as 

a = c a p 2 x2* Vfex P (-r a /r) t (4) 

P = c p Jt, (5) 

Y = C r pX^ 7 Vfexp(-r y /r), (6) 

( 7 ) 

* = ^ Won + ) Vf exp(-r z /n. (8) 


where T is flame temperature, T a , T r T x are activation temperatures, p 
is mixture density, X c is fuel mole fraction, and X 0H and X 0 are mole 
fractions for OH and 0 2t respectively. Numerically, any combinations 
of rate constants (C^ Cp, C r C$ C x ), activation temperatures ( T ^ T r 
T x ), and exponents ( m ^ m Y ) can be adjusted to produce a optimal fit 
between model predictions and experimental data. 

Another model is a one-equation model, based on volume fraction 
only, developed by Khan et al. (1971; 1974). This model characterizes 
soot formation by an Arrhenius -type equation, and the corresponding 
source term is 



where C k is the rate constant, P fu denotes the partial pressure of 
unbumed fuel, <p is the local fuel/air equivalence ratio, and E s is the 
activation energy. The soot oxidation model of Lee et al. (1962) is 
adopted, which is similar to Eq. (9), with ^(p/pj) replaced by f v /^Jf . 
Although the model by Khan et al. is simpler, it seems less physically 
sound than that by Moss et al., since it does not include such 
mechanisms as coagulation which causes decreasing number density 
under constant volume fraction. 


SOLUTION FOR FLAME RADIATION HEAT TRANSFER 

Although it plays an important role, radiation transfer has been 
neglected in detailed combustion analyses due to such difficulties as the 
computational demands and accuracy of various solution methods (Tong 
and Skocypec, 1992), the coupling of radiation and flame structure 
solvers, the accuracy of spectral radiative properties, and efficient 
methods for integrating fluxes and other quantities over the spectral 
range. 

We chose the recently developed YIX method (Tan et al., 1990; 
Hsu et al., 1993) for calculating the radiative heat flux q r and its 
divergence. A solution derived by Hsu and Ku (1994) for finite 
axisymmetric cylindrical enclosures is used for radiation calculations 
from diffusion flames. The YIX method is a numerical approach for 
solving the integral formulation of the radiative transfer equation by 
reducing the order of the multiple distance-angular integrals. The name 
YIX comes from the shape of the pattern of integration points for three, 
two, and four angular directions. One important attribute of this method 
is that these integration points can be pre -calculated and stored, 
significantly reducing computational time. For multi-dimensional 


geometries, discrete ordinate sets are used for angular quadratures. 
Although computationally quite intensive, the YIX method has proven to 
be very accurate and suitable for nonhomogeneous media (Hsu and Ku 
1994; Hsu and Ku 1995). Simplified models such as a fixed percentage 
of local heat loss from equilibrium everywhere, a temperature modified 
from that under equilibrium as T =T eq [\- b{T eq /T eqmax )*] (Leung et 
al., 1991), or a heat sink term given as S $ =-eo{T A -Tl)/l (Kent and 

Honnery, 1987) has each been tested and found inaccurate. We also 
derived a spherical harmonics (P#) approximate solution for nonhomo- 
geneous media following the work of Menguc and Viskanta (1986), but 
found the resulting partial differential equations numerically unstable for 
enclosures involving optically thin regions such as in diffusion flames. 

The radiation solver is coupled to the flame structure solver through 
the energy equation. For diffusion flames, the energy equation takes the 
same form as Eq. (1 ), with <p = H (total enthalpy) and S^=-Vq r , and 

the latter is calculated from a solution for the radiative transfer equation. 
These two solvers, both depending on the temperature, are numerically 
incompatible. The solver for flame structure, soot transport, and the 
energy equation solves parabolic differential equations of Eq. (1) type. 
The solver for radiative transfer, YIX or otherwise, is inherently not 
parabolic, since radiation is an integral phenomenon over ail distances 
and solid angles. We chose an iterative approach and used temperature 
as the convergence criterion. The iteration starts with the flame structure 
solver using an initial guess based on 5, = -£a(7' 4 -£)// to calculate 

velocities, gas and soot concentrations, density, and flame temperature. 
The resulting temperature is used for calculating V q r from the YIX 
solver. The resulting V q r is substituted back into the structure solver 
to update (by averaging the two latest runs) velocities, concentrations, 
density, and temperature. This process is repeated until a convergence 
on the temperature is accomplished between two subsequent iterations 
( AT/ T <0.002). Figure 1 shows a flow chart for the computational 
procedure. For flames calculated here, it typically takes less than 10 
iterations to reach the convergence (Tong, 1995; Ku et al., 1995b). The 
turbulence-radiation interactions, which are shown to be significant by 
Gore et al. (1992), will be included in the future. 

To evaluate spectral radiative properties, the absorption coefficient 
for soot aggregate is calculated from the Rayleigh solution (Ku and 
Shim, 1991) using the refractive index calculated from Habib and 
Vervisch’s (1988) Drude-Lorentz dispersion model. Scattering from soot 
particles is neglected at present to reduce computational time, but will be 
included in the future. The exponential- wide -band model of Edwards 
(1976) is used to calculate the absorption coefficient for C0 2 and H 2 0 
gases. The most accurate way to evaluate the total (i.e„ over the spectral 
range) fluxes and flux divergence is to first calculate them for each 
spectral increment and then calculate the integrals. This is 
computationally inefficient because of the iterations required between the 
flame structure and radiation solver. To focus on the iterations, we 
choose to test the Planck ( a p ) and the Rosseland ( a R ) mean absorption 
coefficients given as (Siegel and Howell, 1992) 


01 [Tl P) = 1TF L av {X ’ T ■ P) e> * {X ’ T) dX ■ 


n de b (T) 


( 10 ) 
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Figure 1. A flow chart for computational procedures of solving flame structure and radiation through iteration. 


RESULTS AND DISCUSSION 

We first compared the radiation heat transfer results from a true 
spectra] integration (1 to 20 fim) against those from mean coefficients 
using pre -calculated temperature and particle and gas concentrations 
distributions for a turbulent ethylene diffusion flame. Figure 2 shows a 
comparison of the soot/gas mixture spectral absorption coefficient at a 
specified condition against both means. Figure 3 shows that the flux 
divergence based on the Rosseland mean and on integration agree quite 
well. We thus base all subsequent calculations on the Rosseland mean 
coefficient. However, it should be noted that the accuracy of radiative 
property values, especially those for soot, requires further study. 

We next tried to identify whether the soot formation model by 
Moss et al. or that by Khan et al. is more accurate by fitting predicted 
flame temperature and soot volume fraction against experimental data for 
a 3.85 cc/s co-flow laminar ethylene flame (Santoro et al.; 1983, 1987). 
The results are shown in Figure 4. The agreements are reasonably good, 
indicating that the models seem to be physically sound. Values for rate 
constants in Moss' model are listed in Table 1. For Khan’s model, we 
used C k = 16.8 kg/Nms, m k = 3, and E s = 40000 cal/mole. Although 
both models yield similar fit to experimental data, other calculations 
reveal that Khan's model is overly sensitive to temperature and tend to 
over-predict centerline volume fraction distributions. We therefore chose 
Moss' model for all subsequent calculations. Figure 4(a) also shows that 


the inclusion of radiative transfer does provide a more accurate 
magnitude and trend in temperature profiles than adiabatic values. 
Figure 4(c) shows that the temperature converges smoothly in less than 
10 iterations. Both models yield fairly accurate predictions for gas 
velocity (Tong, 1995). It should be pointed out that, except for the work 
by Fairweather et al. (1992), other works on soot modeling (Moss et al., 
1988; Syed et al., 1990; Kennedy et al., 1990; Kennedy et al., 1991) 
made comparisons against soot volume fraction data only. We base our 
comparisons on ail available gas velocity, temperature, and soot volume 
fraction data. 

We then examined the generality of Moss’ model by fitting 
predictions to experimental data. We used two groups of data, one from 
Santoro et al., (1983. 1987) with the same fuel and burner but at different 
flow rates, and the other from Greenberg and Ku (1996b) for 1-g and 0- 
g 2.30 cc/s laminar jet acetylene names. Table 1 lists burner dimensions 
and flow conditions for these flames. Table 2 lists values of rate 
constants which produce a reasonable fit for soot volume fraction 
between model predictions and data, as shown in Figures 5 through 10 
We adjusted all rate constants (C* C r C* CJ to find a reasonable 
fit. Even though activation temperatures {T a , T r TJ and exponents 
( m or m y) can a lso be adjusted for possibly better fits, we chose to use the 
same values as those in the original works, since our experience indicates 
that these parameters are less influential (Tong, 1995). 
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Figure 2. Absorption coefficients for soot/gases mixture at 1400K, with 
0.1 mole fraction for CO 2 and H^O, and a soot volume fraction of to 6 . 



Figure 3. The flux divergence contours for a simulated ethylene jet flame 
(D = 0.58 mm, Q = 3.96 cm 3 /scc, and Re = 536). 



265 

IV-5 





Table 1 Burner Dimension and Flow Conditions for Flames in This Study 


Fuel 

Dfuel 

(mm) 

D*r 

(mm) 

Q-fuel 

(cc/s) 

V 

(cm/s) 

^ e D/utl 

Qcr 

(cc/s) 

(cm/s) 

Ethylene (C 2 // 4 ) 

11.1 

101.6 

2.30 

2.38 

16.3 

713.3 

8.90 

Ethylene ( C 2 H A ) 

11.1 

101.6 

3.85 

3.98 

27.2 

713.3 

8.90 

Ethylene (C 2 H A ) 

11.1 

101.6 

4.90 

5.06 

34.7 

1068.3 

13.3 

Acetylene {C 2 H 2 ) 

1.65 

NA 

2.30 

107.7 

179.4 

NA 

NA 

Ethylene (C 2 H A ) 

3.0 

NA 

36757 

5200 

11725 

NA 

NA 

Acetylene ( C 2 H 2 ) 

0.508 

NA 

10.0 

4933.8 

2550 

NA 

NA 


Table 2 Values of Soot Formation Rate Constants for Flames in This Study 


Flames 

C a x 10- 7 

C Q X 10^ 

C r x 10 1S 

Q > 

C x x 10‘ 2 

2.30 cc/s C 2 H A (1-g) 

0.7 

1.0 

0.62 

144 

1.4 

3.85 cc/s C 2 H A (1-g) 

2.0 

1.0 

0.62 

144 

1.4 

4.90 cc/s C 2 H 4 (1-g) 

3.0 

1.0 

0.62 

144 

1.4 

2.30 cc/s C 2 H 2 (1-g) 

5.8 

1.0 

0.62 

144 

3.62 

2.30 cc/s C 2 H 2 (0-g) 

2.7 

1.0 

0.62 

144 

2.17 

2.30 cc/s C 2 H 2 (1-g) 

[m a = my- 4] 

320.0 

1.0 

39.7 

144 

2.41 

2.30 cc/s C 2 H 2 (0-g) 

(m a = rrty= 4] 

230.4 

1.0 

39.7 

144 

2.16 

0.61 1pm C 2 H A (1-g; 

0.2275 

1.3 X 10 3 

2.0 x 10 13 

5040 

0.002 

10.0 cc/s C 2 H 2 (1-g) 

0.2275 

1.3 X 10 3 

2.0 x 10 13 

5040 

0.002 


Note: SI units. In aLl cases, T a = 46100 fC T y - 12600 K, and T % - 
19680 K. Unless specified otherwise, m a - m y = 1. 


Numerically, we feel the fits can be improved, but we would like to 
examine the aspect of improving the models themselves first. As can be 
seen from the figures, the 4.90 cc/s ethylene and both 1-g and 0-g 2.30 
cc/s acetylene flames are smoking (i. e., releasing soot particles from the 
flame tip), while others are non-smoking. In Figures 5 and 6, for 4.90 
and 2.30 cc/s respectively, we included a dashed curve for predictions 
using the same rate constants as those for 3.85 cc/s. This shows that 
Moss model is not sensitive enough to changes in fuel flow rate. 
However, this can be easily improved by adjusting C ^ the nucleation 

rate constant. The fact that the fit looks better for 2.30 cc/s flame (Figure 
6) than for 4.90 cc/s flame (Figure 5) suggests that further adjustments 
of rate constants may be required going from non-smoking (2.30 and 
3.85 cc/s) to smoking flames (4.90 cc/s). 

When applied the models to 1-g and 0-g laminar jet acetylene 
flames, we first tried to fit the data by adjusting the rate constants C a and 


C x only, while keeping m a = m y = 1. The results are shown in Figures 
7 and 8. The fits are obviously not as good as those for co-flow ethylene 
flames, and perhaps expectedly so due to the fundamental differences 
between jet flames and co-flow flames. However, further adjustments 
on model coefficients, including setting m a - m y - 4, yield better fits, as 

shown in Figures 9 and 10. Notice that the fits to 1-g and 0-g flame are 
obtained with only minor adjustments to C a and C r 

Based on the facts that the model can produce reasonably accurate 
fits to data for flames of different types, flow rates, fuels, and l-g/0-g 
conditions, we feel that Moss’ model is physically sound. Some 
patterns about the effects of each model coefficient on the soot volume 
fraction distribution are identified as follows. 

1. The overall shape and peak locations of volume fraction distribution 
are fairly independent of all rate constants (C^ CfiCrCs.C z ). 

2. An increase in nucleation ( C ^ or surface growth (C r ) rate constant 

will cause an overall increase in soot volume fraction, with higher 
sensitivity to C^than to C a . 

3. Coagulation (Cg) and oxidation (C z ) rate constants have negligible 
effects on volume fraction distributions near the burner exit, but 
have significant effects at higher locations in the flame. An increase 
in C x will cause a decrease around the edge, and an increase m Cp 

will cause a decrease around the center region of the flame. 

4. An increase in both exponents ( m a , m y ) will cause an increase in 
soot volume fraction around the centerline region at higher locations 
in the flame, and will shift the soot volume fraction peak locations 
inward (toward fuel-rich side, as if the soot yield were dependent on 
an intermediate species rather than the parent fuel). 

5. An increase in thermophore tic velocity v ; will cause an increase in 
soot volume fraction around the centerline region. 

Finally, we tested the model on turbulent jet flames. We compared 
against the published data of Kent and Honnery (1987) for a 1-g ethylene 
flame, and against our data for a 1-g acetylene flame. Figure 11 shows 
comparisons of two model predictions, one using the YTX solver and the 
other using the simple sink term, against Kent and Honnery's (1987) 
data. Figure 12 shows predicted volume fraction maps for 1-g and 0-g 
acetylene flames, and a comparison against our 1-g data. Again, the 
model was able to predict the correct level and overall shape of volume 
fraction and temperature distributions, even though we only made 
limited adjustments and used the same rate coefficients for both flames, 
due to the excessive computational time required and modeling 
uncertainties in turbulence -radiation interactions and radiative properties. 

In Figure 11, the discrepancy may be caused by the missing turbulence- 
radiation effects. In Figure 12, we can modify- the oxidation constant 
(C x ) to get a better fit. The burner dimension and flow conditions and 

rate coefficients for these flames are listed in Tables 1 and 2. The 
noticeable differences in rate coefficients between turbulent and laminar 
flames are consistent with Moss’ published values. 

We plan to adopt more accurate means of evaluating total radiative 
quantities, such as the wide-band model used by Song and Viskanta 
(1987). Their derivation of a ’’mean" radiative transfer equation to include 
turbulence-radiation effects is also worth noting. A joint pdf (Gore et al., 
1992; Janika and Kollmann, 1987)) could be developed to correlate 
enthalpy, temperature, and mixture fraction for modeling turbulence- 
radiation effects,. Turbulence-buoyancy interactions (e. g., Fairweather et 
al., 1992) will be examined for 0-g flames. 
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Figure 12. Comparisons of predictions to data for a 1-g 10.0 cc/s turbulent acetylene jet diffusion flame. 
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